Ergodicity and Slowing Down in Glass-Forming Systems with Soft Potentials: No 

Finite-Temperature Singularities 
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The aim of this paper is to discuss some basic notions regarding generic glass forming systems 
composed of particles interacting via soft potentials. Excluding explicitly hard-core interaction we 
discuss the so called 'glass transition' in which super-cooled amorphous state is formed, accompa- 
nied with a spectacular slowing down of relaxation to equilibrium, when the temperature is changed 
over a relatively small interval. Using the classical example of a 50-50 binary liquid of -/V particles 
with different interaction length-scales we show that (i) the system remains ergodic at all temper- 
atures, (ii) the number of topologically distinct configurations can be computed, is temperature 
independent, and is exponential in N. (iii) Any two configurations in phase space can be connected 
using elementary moves whose number is polynomially bounded in AT, showing that the graph of 
configurations has the 'small world' property, (iv) The entropy of the system can be estimated at 
any temperature (or energy), and there is no Kauzmann crisis at any positive temperature, (v) The 
mechanism for the super-Arrhenius temperature dependence of the relaxation time is explained, 
connecting it to an entropic squeeze at the glass transition, (vi) There is no Vogel-Fulcher crisis at 
any finite temperature T > 0. 
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I. INTRODUCTION 

It is not uncommon to read in papers devoted to the 
glass-transition a statement of the type 'it is well known 
that glass forming systems lose ergodicity'. It is even 
more common [l[ to use fits of the a relaxation-time r a 
in such systems to the Tamam-Vogcl-Fulcher formula 

(1) 

which indicates a belief that the relaxation time actu- 
ally diverges at some temperature Tvf > 0. Related to 
these issues is the concept of the Kauzmann temperature 
which is a finite temperature Tk > where the ex- 
trapolated entropy appears to vanish. In this paper we 
argue that these related questions should be addressed 
with care; we wish to clarify which of these issues can 
actually appear in glass-forming systems, and which of 
them are only a consequence of inadequate simulations 
or interpretations, or even of a confusion of questions. 

In some parts this paper is a review, or an outgrowth 
of our own efforts to understand basic concepts of glass- 
formation, which became more transparent in some of 
our earlier work [H, [J, 5], and whose findings we combine 
in this paper. Glasses and their formation have occupied 
researchers for many decades, with ideas being first devel- 
oped on theoretical bases @, and, in more recent years, 
being studied extensively with the help of computer sim- 
ulations 0, H, M, E3- Such simulations vary from very 
realistic models to toy models with simplified dynamics, 
but correspondingly faster simulations. We can summa- 
rize our findings in the following logical structure: 

1) One has to distinguish between large (but finite) 
systems as compared to infinite systems. Infinite systems 
pose difficult conceptual problems, because in this case, 
density and close packing are not tightly related (any 
small lowering of density will allow for arbitrarily large 
voids, and the jamming problems disappear). 



We will therefore focus our discussion only on systems 
with a finite, but arbitrary, number of particles. In this 
view an Avogadro number of particles is finite, and it is 
important to state this fact. 

2) There is an essential difference between systems or 
particles with a hard core and systems of particles inter- 
acting with soft repulsive potentials. In the first case, 
there is obviously a density where the particles can not 
move any more (or at least not all of them can be moved) . 
This poses interesting problems of jamming, or contact 
geometry. These issues have been studied in depth by 
Stillingcr and his school, with a careful analysis of differ- 
ent types of 'movability' We point out that these 
problems, while very interesting from the geometrical 
point of view (see also the book by Conway and Sloanc 

do not really address the questions of what a generic 
glass is. We will thus focus attention to systems with a 
finite number of particles, and with soft potentials. 

3) There is another issue which deserves attention: 
time. Many concepts of glasses are based on the notion 
that if something does not happen before a given time, 
then it will never happen. Such reasoning is unaccept- 
able from a theoretical point of view, since a natural time 
scale does not exist. We will argue that in systems with 
a finite number of particles interacting via soft potentials 
there is no singularity in the relaxation time at any pos- 
itive temperature. In particular such systems undergo a 
glass 'transition' in the sense that their relaxation time 
increases without limit as temperatures are lowered, but 
they remain ergodic. Observing the consequences of er- 
godicity may necessitate waiting for an unbounded, but 
still finite time. At any temperature larger than T = all 
of phase space is available to the system, all the config- 
urations are dynamically connected; the configurational 
entropy can be computed, yielding a finite result for any 
given finite energy (or finite temperature). We stress 
and reiterate that this paper deals with systems having a 
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finite number of particles, which are interacting via soft 
potentials and are being observed for an unbounded time. 

One can exemplify the discussion with the help of many 
simple models, and for concreteness we choose the classi- 
cal example of a glass forming binary mixture of particles 
interacting via a soft 1/r 12 repulsion with a 'diameter' 
ratio of 1.4. More or less related models can be found 
in [H, EI EE EE El El- This 2-dimensional model 
had been selected for simulation speed and, more impor- 
tantly for the ease of interpretation. We refer the reader 
to the extensive work done on this system [E 0, 0, [E El • 
It shows that it is a bona fide glass-forming liquid meeting 
all the criteria of a glass transition. In short, the system 
consists of an cquimolar mixture of two types of parti- 
cles with diameter (T2 = 1.4 and o\ = 1, respectively, 
but with the same mass m. The three pairwise additive 
interactions are given by the purely repulsive soft-core 
potentials 



u a b = e ( — ) , a, 6=1,2 



(2) 



where a aa = a a and a a b = {o~a + ci,)/2. The cutoff radii 
of the interaction are set at A.5o~ a b- The units of mass, 
length, time and temperature are m, eri, r = aiy/m/e 
and T = e/ks, respectively, with ks being Boltzmann's 
constant. 

It has been shown that for temperatures T > 0.5 the 
system behaves like a liquid. For temperatures lower 
than 0.5 the system begins to slow down, the correlation 
functions stop decaying exponentially, and can be fitted 
within reason to a stretched exponential form, cf. Fig. [T] 
The time constant (or so-called a relaxation time) can 
be fitted to the Vogel-Fulcher form fl} for temperatures 
not too close to Tvf ■ The model shows the expected be- 
havior of the specific heat at the temperature range that 
is considered the 'glass transition' (i.e., where the relax- 
ation time increases rapidly), and the 'super- Arrhcnius' 
dependence of the viscosity on the temperature. In short, 
the system appears as a good example to consider to elu- 
cidate the fundamental issues that concern us here. 

In Sect. [II] we explain how to discretize the configura- 
tional space using Voronoi tesselations and their dual De- 
launay triangulations. This offers the basis of the demon- 
stration of ergodicity at any finite temperature. We also 
review the mathematical results concerning the size of 
this phase space (exponential in the number of particles 
N) and the fact that its graph has the small-world prop- 
erty in the sense that any configuration can be reached 
from any other using a polynomial number of steps. In 
Sect. IIIII we discuss the configurational entropy, and ex- 
plain that there is never a Kauzmann crisis in this system. 
Finally in Sect. lIVl we present the fundamental reason for 
the ' super- Arrhenius' dependence of the relaxation time 
on the temperature. This is due to the fast decrease in 
the concentration of some quasi-species, leading to an en- 
tropic squeeze. In Sect. [V] we summarize the paper and 
offer concluding remarks. 
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FIG. 1: (Color online:). Typical correlation functions at dif- 
ferent temperatures for the model under discussion. Shown 
are the shear stress-shear stress correlation function which de- 
cays exponentially at high temperatures, and then, at about 
T — 0.5, changes character to a stretched exponential form. 
At temperatures lower than 0.4 the correlations decay to a 
finite value, which is determined by the mean shear modulus, 
indicating that the material behaves like a solid for the given 
measurement time. For further details about these measure- 
ments see 12011. 



II. VORONOI CONSTRUCTION AND 
ERGODICITY 

A. Voronoi tessellation and Delaunay triangulation 

We begin by discussing the possible configurations of 
the system in a systematic way, using the time-honored 
Voronoi polygon construction [7fl. It associates with ev- 
ery configuration of the particles a subdivision of position 
space into polygons, one per particle. These polygons 
will also be called cages. The polygon associated with 
any particle contains all points closer to that particle 
than to any other particle. The edges of such a poly- 
gon arc the perpendicular bisectors of the vectors joining 
the center of the particle (actually the coordinates of the 
point particles we consider). As had been noted in 
when periodic boundary conditions are used, the average 
number of sides of the polygons is exactly 6. This fol- 
lows because the Euler characteristic on the torus is 0: 
X = V — E + F = 0, where V, E, and F are respectively 
the numbers of vertices (corners), edges and faces in the 
given polyhedron. 

Typical such Voronoi tessellations for two different 
temperatures are shown in Fig. O Since there are two 
types of particles (small or large, or blue and red (respec- 
tively)), in order to have a mapping between the parti- 
cle positions and the number of sides in the polygons of 
the Voronoi tessellation, we distinguish between polygons 
having a small and large particles in their center. Thus, 
a coloring scheme of cells will take into account not only 
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FIG. 2: Upper panel: The Voronoi polygon construction in 
the liquid state at T = 3 with the seven-color code used 
in this paper. Small particles in pentagons (heptagons) are 
light green (dark green) and large particles in pentagons (hep- 
tagons) are violet (pink) . Lower panel: a typical Voronoi con- 
struction in the glass phase at T = 0.1. Note the total dis- 
appearance of liquid-like defects (large particles in pentagons 
and small particles in heptagons). 



the number of sides, but also the type of the particle (big 
or small) in each cell. The tessellation obtained without 
distinguishing the type of particles will be referred to as 
'colorless' and so, below, will be the resulting triangu- 
lation. When distinguishing the two types of particles, 
the tessellation and its triangulation will be referred to 
as colored. 

Note that in a generic configuration the Voronoi de- 
composition has the property that exactly three lines 
meet at a point. Accordingly, the dual graph of such a de- 
composition is a triangulation. Such triangulations are 
called Delaunay triangulations. They are obtained by 
simply connecting the centers of those particles which are 
separated by one edge of the polygon. In other words, 
the vertices of the triangles coincide with the particles 
at the centers of the Voronoi cells. Let us remark that 
this construction can be also done in 3-dimensions (or 
more), where now instead of edges of polygons we have 
faces of polyhcdra and four faces will meet at a point. 
In this case, the dual graphs are now made of tetrahc- 
dra. In 3-dimcnsions we will call the graphs obtained 



in this way Delaunay tetrahedrization. It is noteworthy 
that there is an important difference between 2 and 3- 
dimension. In 2-dimensions, the points of any triangula- 
tion can be arranged in such a way in position space that 
the triangulation is indeed a Delaunay triangulation, as 
described above, that is, the dual of a Voronoi decom- 
position. However, in 3-d there exist tetrahedrizations 
which arc not duals of Voronoi tessellations, or in other 
words, cannot be obtained starting with particles spread 
in 3-dimensions, ma ppe d onto a Voronoi tessellation and 
then tetrahedrized ([2l|). 

Having outlined the description of configurations in 
this model (and in fact all models of this general type), 
we consider now the question of dynamic accessibility, 
viz. ergodicity. The model under consideration is Hamil- 
tonian, but in simulations it is coupled to a heat bath 
of temperature T . There are several methods for achiev- 
ing such a coupling, leading to the study of this model 
in either N . V, T or N , p, T ensembles. In either en- 
semble the mean energy per particle is fixed, but there 
can be arbitrarily large fluctuations in the energy. It is 
precisely these fluctuations which are the source of ergod- 
icity. The energy is partially kinetic and partially poten- 
tial. By ergodicity, one means that the natural motion of 
the system can reach every point in phase space. For sys- 
tems having one type of (indistinguishable) particles, it is 
not important to look at particles exchanging positions: 
indeed every point in phase space can be obtained by 
moving particles within their cages and translating cages 
as is necessary. In contrast, when there is more than 
one type of particle (here we have two types) , ergodicity 
means that particles must also be able to change their 
relative positions in all imaginable ways. It is therefore 
useful in both cases to distinguish between thermal mo- 
tion of individual particles rattling in their own 'cages', 
without really changing their relative positions, and large 
scale movements and exchanges of particles that change 
the macroscopic configuration. Only the second class of 
movements is important in terms of the configurational 
entropy of the system, which arises from counting dif- 
ferent configurations, for which the thermal motion does 
not matter. For example, in systems with hard core in- 
teractions, there are various levels at which motion can 
be restricted, as we have mentioned before. But here, we 
consider soft particles. 

Since we are starting from a problem with particles 
which live in a concrete space, we can focus on the re- 
stricted class of Delaunay graphs (which actually coin- 
cides in 2-d with all triangulations, as we have explained 
earlier). In other words, in 2-dimensions (3-dimensions) 
the natural phase-space for our system is the set of all 
Delaunay triangulations (tetrahedrization). 

In 2-dimensions the elementary change in the Voronoi 
tessellation is obtained by a so called T\ process as seen 
in Fig. [3] (This operation has various names in the 
literature: Gross- Varstedt move, Pachner move, flip.) 
[H [H, El EE El- Note that thc operation of a T Y 
process in the Voronoi tessellation translates to a flip in 



FIG. 3: A Tl process in the Voronoi tessellation (in red lines) 
and the equivalent flip process in the Delaunay triangulation 
(dashed curves) 



the Delaunay triangulation, see Fig. [3] 

To simplify things further, we consider instead of the 
triangulations of the torus the triangulations T of the 
sphere. The simplification is that the genus is 0, and that 
more is known about the combinatorics of triangulations 
of the sphere than of the torus. Given N, correspond- 
ing to the number of particles in the original model, we 
let T denote a triangulation of the sphere with N nodes 
and we let IV, o denote the set of all such triangulations. 
By this, one means the set of all combinatorially distinct 
rooted simplicial 3-polytopes. In particular, a triangu- 
lation should not have any 'double edges'. We further 
refine the definition, by distinguishing 2 types of nodes 
in the triangulation: We first number the nodes from 1 
to N and then define 2 types of nodes. Those with even 
index are the 'small particles' and those with odd index 
the 'large' ones. This means that the triangulation has 
the same number of large and small nodes (up to a dif- 
ference of one). We will also call the two types of nodes 
two colors. 

We shall call the odd nodes blue and the even ones 
red and will refer to the triangulations as colored trian- 
gulations. Once the colors arc assigned, the numbers are 
again forgotten. The set of all colored triangulations with 
N nodes will be called IV ■ This is our phase space and 
the dynamics is mapping points in this phase space to 
other points. Having understood that, we can now ask 
the major question of this section: 

The question of ergodicity: is every point in phase- 
space accessible through a sequence of motions that are 
flips. The answer is not obvious from the outset. For ex- 
ample in 2-dimensions one can consider a triangulation 
that includes locally a form as in Fig. 0] Obviously, a 
form like this is stuck, since one cannot do any flip with- 
out doubling one of the links. (One can see this easily 
by noticing that the configuration in Fig. [4] is a valid 
triangulation of the sphere, and it is the unique trian- 
gulation with four points.) The question of ergodicity is 
even less trivial in 3-dimension, where there is a whole 
class of unflippable tetrahedrizations. However, they are 
not duals of Voronois. But, the relevant question (in any 
dimension) is whether the triangulations (tetrahedriza- 
tions,. . . ) which are duals of Voronoi tessellations are 



FIG. 4: A configuration with four nodes that cannot be 
flipped, because that would lead to double links. In 3- 
dimensions there is more than one corresponding graph, but 
none of them is the dual of a Voronoi [2ll [. 



indeed all connected in one large component. The an- 
swer is nevertheless 'yes' in both 2 and 3-dimcnsions as 
we explain next. 



B. Demonstration of ergodicity 

One can attack this problem in two ways, the first be- 
ing more physical, and the second more mathematical. 
The physical argument is trivial. Since the particles are 
soft, and the energy is not bounded, they can be moved 
around each other in any way one can imagine. We stress 
again that while this needs perhaps a lot of energy, the 
large fluctuations of the energy will guarantee that this 
will eventually happen, rarely, but surely. And hence er- 
godicity is obvious. The only care one must observe is 
that moves through degenerate situations (4 lines meet- 
ing at a Voronoi vertex) must be avoided. 

The question is more intriguing when formulated in 
terms of triangulations alone, since we have seen in Fig. 
2] that a local configuration in which 3 links emanate from 
a node cannot be flipped. Could it be that the triangu- 
lation is so imbricated that in fact none of the edges can 
be flipped? Indeed, for 4 particles this is exactly what 
happens, but then, the phase space consists of exactly 
one triangulation and thus there is no need to move any 
link. In 1936, Wagner (2?| showed that a finite number of 
flips will transform any triangulation of the sphere into 
a 'Christmas tree,' which is the configuration shown in 
Fig. [5] See also [28| for a discussion of several related 
issues. Since one can undo the flips, this implies immedi- 
ately that any two triangulations can be connected by a 
sequence of flips going through the Christmas tree. There 
is abundant literature on this question [2l|, [24[ which also 
plays a certain role in the classification of 3-manifolds. 
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asymptotically 



FIG. 5: The Christmas tree. Two nodes are at the bottom, 
the others (only 4 shown) are in the stem of the tree. 



C. The size of phase space 

The possible states of our system of triangulations with 
N nodes is the set Tjv of all possible colored triangula- 
tions. The set Tjv has, as we will see, a number of el- 
ements which grows like C N for some constant C > 1. 
It is thus a discrete space with a finite number of states. 
To describe the dynamics of flipping in a geometric way, 
one should view this set as the dynamical graph Q, whose 
nodes are now the elements of the set Tn (not to be con- 
fused with the nodes (particles) of any triangulation T) 
and two of its nodes are linked if one can be reached 
from the other by a flip. (This makes an undirected 
graph, since one can flip back and forth.) The reader 
should note that there are two graphs in this discussion: 
Each triangulation is a graph with N nodes, and 37V — 6 
links (by Euler's theorem), while the graph Q has about 
C N nodes, and about 3N — 6 links per node. This last 
statement follows because in every state of Tjv, one can 
choose which of the 3N — 6 links of the triangulation T 
one wants to flip. However, there will, in general, be 
somewhat fewer links which are candidates for flipping, 
because whenever there is a node of degree 3 in the tri- 
angulation T its links can not be flipped (a tetrahedron 
is unflippable). 

Finally, given any two elements in Tjv, that is, any two 
triangulations with N nodes, we will show below that 
0(N 2 ) flips are sufficient to walk on the graph Q from 
one to the other. Thus, the diameter of the graph Q is at 
most 0(N 2 ) while it has 0(1) N vertices. This means that 
Q has the 'small- world' property (29[. It has also small 
clustering coefficient, since there are very few triangles 
in the graph Q (it is difficult to get from a triangulation 
back to the same triangulation with 3 flips). 

In the remainder of this section, we spell out these 
statements. They are well-known for uncolored graphs, 
so the only task is to prove them for the colored graphs, 
see @. 

We first state two known results for the set T/v,o of 
uncolored triangulations: 

Lemma 1 23l . 1241 . |30[ The number of elements in IV o is 



256\ JV " 3 3 



27 



16V67riV 7 



(3) 



The distance between any two uncolored triangulations 
is at most 6N — 30 flips. 

For the case of the colored graphs, with iV rec j = iVbi U e + 
c and c € {0, 1}, that is, about equal number of red and 
blue nodes, one has 

Lemma 2 The number of elements in Tjy is asymptoti- 
cally bounded above by 



N-3 



2 n (256 

V 27 7 16V67riV 5 



(4) 



and below by the expression ([3|). The distance between 
any two colored triangulations in Tjv is bounded by 



dN 2 + C 2 



(5) 



flips with some universal constants C\, C%. 

We note that the phase space as defined here is (obvi- 
ously) independent of the temperature. We can thus con- 
clude that the present classical model of glass-formation 
does not suffer from any issue of loss of ergodicity. Ac- 
cordingly, it should have a valid statistical mechanics at 
any temperature T > 0. Next we show that indeed its 
configurational entropy S never suffers any finite tem- 
perature crisis, and there is no Kauzmann temperature 
where S ^ 



III. CONFIGURATIONAL ENTROPY 
A. Statistical Mechanics 

Needless to say, the configurations discussed in the pre- 
vious section have different energies and therefore the 
probability to see any particular one can be strongly de- 
pendent on the temperature. To discuss the temperature 
dependence of the configurational entropy of this system 
we need to review the statistical mechanics that was in- 
troduced for this system in Q]. The basis of the analysis 
is again the Voronoi tessellation. 

To get some systematics, we first define a 'typical' en- 
ergy for each type of topological cell, with an average 
taken over all cells with a given number of sides and a 
given particle type (big or small) in its center. The aver- 
e ( \ 12 

age is then e; = ( Ylk=i 6 ( ) ) wnere Ei is the num- 
ber of edges associated with that cell type (i), and r,fe 
being the distance to the particle in the adjacent Voronoi 
cell and the average is over all particles of the same type i. 
In Fig. [S] we present the values of these energies measured 
numerically as a function of the temperature, following 
a protocol of slow cooling. In the range of temperatures 
explored by simulations there are 10 different cell types, 
(large particle or small particle in squares, pentagons, 
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FIG. 6: (Color online). The average energies of the Voronoi 
cells as a function of the temperature as measured in the 
simulations. 



hexagons, heptagons and octagons), but octagons and 
squares are not shown since they already disappear at 
relatively high temperature, much above the glass tran- 
sition. 

We learn from these data graphs that the different cell 
types have clearly split energies throughout the interest- 
ing temperature range, and that these energies are only 
weakly dependent on the temperature. Within the tem- 
perature range of interest we can focus on the six types of 
cells; denote by {Ni}® =1 the number of cells of each type, 
with number of edges E^, ordering them by the mean 
energy ej with i = 1 being the highest (large particle in 
a pentagon) and i = 6 being the lowest (small particle 
in an heptagon). Additional important properties of the 
cell types are their areas fij and their shapes; the first 
affects the enthalpy term and both affect the configu- 
rational entropy when we count the number of possible 
tilings of the plane. 

With this in mind we can construct the statistical me- 
chanics of this system by considering the free energy 
G = U +pV—TS. We should stress at this point that one 
could aspire for more accurate statistical mechanics, con- 
sidering for example not only the type of particle inside 
the Voronoi n-gon, but also who are the neighbors (small 
or large particles) (which resembles somehow the plaque- 
tte expansions in statistical mechanics). Such a choice 
would have allowed a better treatment of the tendency 
of hexagons with large (or small) particles for example 
to crowd together to minimize the pV term in the free 
energy. The price is that the number of quasi-species 
increases to 42 (6, 7 and 8 for each pentagon, hexagon 
and heptagon respectively) . While doable and a bit more 
precise, this more involved statistical mechanics does not 
shed more light on the issues of principle that interest 
us in the present paper, and therefore we do not discuss 
such improvements any further. 

Coming back to the simpler variant, we note that in 
the free energy the value of U is J2^=i^i e i- The pV 
term is simply p^] i=1 iVjfij. Lastly, we need to estimate 



FIG. 7: The concentrations of all the Voronoi cells as a func- 
tion of T. Note the close-to-degeneracy of pairs of cells: the 
upper pair are the glass-like quasi-species, the middle pair 
are the two hexagons and the lower pair the liquid-like quasi- 
species. Observe the existence of two typical temperatures 
ranges, one around T2 where the liquid-like quasi-species de- 
plete severely, and the other around Ti where the hexagons 
deplete rapidly. 



the entropy term. In principle this should be computed 
from the number of possible complete tilings of the area 
by Ni cells of each type with its given area and shape, 
subject to the Euler constraint J2i=i NiEt = 6N, where 
Ei is the number of edges of the i'th polygon. This is a 
formidable problem. A useful estimate can be obtained 
by considering the area only, and filling space starting 
with the largest objects, then the next largest, and so 
on, until the smallest are fit in. We do this by dividing 
the (remaining) volumes into boxes, and studying the 
combinatorial filling of these boxes. Denoting the possi- 
ble number of boxes to fit the largest cells by N\ = V/ili , 
then the number of boxes available for the second largest 
cell by A/2 = (V — NiQi)/Q,2 etc., the number of possible 
configurations W is 



w = n 



(6) 



Using the abbreviation x% = Ni/Mi we compute directly 
xi = CjClj/ Ylj—j CjClj where Cj is the number concentra- 
tion of each defect. We can now compute S = fee In W 
and write G together with a Lagrange multiplier for the 
Euler constraint, 

6 6 6 

G = 12 NiCi + p 12 NiQi + x 12 N * Ei 

i—1 i—1 i—1 

6 

+T^N k [x k hix k + {l-x k )\n(l-x k )] . (7) 
fc=i 

The chemical potential fj,i = dG/dNi is then, for i = 
1,...,6, 



k=l 
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We now recognize that Voronoi cells of different values 
Ei but with the same size particle (small or large) are 
in equilibrium, each one being able to change to an- 
other, but small particles cannot change to large par- 
ticles, and therefore in equilibrium there exist only two 
independent values of fa, one for the small particles (j,s 
and one for the large particles /^l, and we have 9 un- 
knowns - six values of c^, 2 values of n and one Lagrange 
multiplier A. This is precisely balanced by the 6 equa- 
tions ©, the Euler constraint, and the two constraints 

Si=i c i = Sfc=4 c k = 1/2- These equations could be 
solved numerically using the precise values of Oj(T) and 
Ci(T) as measured in the simulation. The approximate 
calculation of the entropy however does not warrant such 
a detailed calculation. In reality, calculating the average 
areas of the cell types in the numerical simulations, wc 
discover that to an excellent approximation these fall in 
two classes, smaller cells of area fig when small particles 
are in them, and larger cells of area Oz, where large par- 
ticles are enclosed. These areas again are only weakly 
dependent on the temperature. Then the whole system 
of equations simplifies to two analytically tractable sets 
of equations 

p, L = a+T In a + XEi , {* = 1, 2, 3} , 

}i s = ei + Thia + XEi , {i = 4,5,6}, (9) 

together with the above mentioned three constraints. In 
p, we have absorbed terms that added to fi in this special 
case. 

In Fig. we show the solutions of these equations when 
we use values of ej taken from the Fig. [5] at T = 0. We 
learn from these results that the statistical mechanics 
predicts that the first quasi-species to disappear are the 
large particles in pentagons and the small particles in 
heptagons. While the first is the highest in energy, the 
disappearance in tandem of the second is a result of the 
Euler constraint, and could not be guessed a priori. In 
previous work the hrst disappearing quasi-species were 
called 'liquid-like', since they are common in the liquid 
state and their concentration is exponentially small in 
the glass state. The region of temperature where their 
concentration falls off rapidly was identified with the re- 
gion of slowing down. In fact, in Q a quantitative rela- 
tion between the concentration of these liquid-like quasi- 
species and the relaxation time was derived, explaining 
the slowing down as a result of an entropic squeeze. We 
will return to this issue in Sect. HVl 

The statistical mechanics predicts a second transition 
(cf. Fig. [7]). Below some temperature the concentration 
of hexagonal cells is predicted to be exponentially small, 
and the system retains only pentagons with small par- 
ticles and heptagons with large particles. These are re- 
ferred to as the 'glass-like' quasi-species. Indeed, in [i[ 
it was found that a phase made of only pentagons with 
small particles and heptagons with large particles exists 
and is stable at low temperatures, see Fig.[5J Upon warm- 
ing up such a phase, at a temperature roughly around 
Ti, a sizable number of hexagons appears to form the 




FIG. 8: (Color online). An example of a stable phase formed 
by glass-like quasi-species at very low temperatures. 



generic glassy state. Upon further warming, crossing 
a temperature roughly around T2, a sizable number of 
liquid-like defects brings the system to a liquid state. The 
actual values of T\ and T2 can be understood from this 
model. Denote by ct, ch and cg the concentrations of 
the liquid-like, hexagons and glass-like quasi-species, and 
by ei = ei + eg as 12.48 as the energy associated with the 
liquid- like quasi-species, by ejj = e 2 + e 5 ~ 11-94 the en- 
ergy of the hexagons, and eg = £3 + 64 ~ 11.76 the energy 
of the glass-like defects. The theory predicts that ratios 
C£/ch and ch/cq are of the order of exp[— (q — ejj)/T\ 
and exp[— (ch — e c)/T] respectively. As an estimate of 
T2 and T\ take these ratios to be, say, of the order of 
1% ~ exp(— 5) and observe that such ratios are obtained 
for T = T 2 as 0.11 and T = T x as 0.04. It is important to 
notice that eh — eq could be positive rather than nega- 
tive, and then the system would crystallize on a hexag- 
onal lattice. Such a lattice can exist in this system only 
when the particles phase separate into two pure hexago- 
nal lattices of small and large particles respectively, with 
an interface in between. Such a phase may even be the 
ground state, but seems to be inaccessible in dynamical 
experiments starting from random organizations of small 
and large particles. 



B. Nonexistence of Kauzmann temperature 

Assuming that indeed the realized state at T = is 
the state shown in Fig. [51 but equivalently if the ground 
state were made of hexagons, it is obvious now that any 
finite temperature will allow the appearance of the other 
quasi-species (hexagons inside the phase in Fig.[8jor pairs 
of pentagons and heptagons in the hexagonal phase if the 
latter were the ground state). Focusing on the first sit- 
uation, we understand that any pair of hexagons costs 
a given amount of energy eh which is a given A above 
the ground-state eg- The calculation given by Eq. (J7]) is 
legitimate at any temperature T, and the configuration 
entropy is approximately correct as stated there. One 
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can improve the calculation of the entropy compared to 
the approximation employed above, but there is nothing 
extraordinary that is expected at any temperature. We 
thus state that the configurational entropy is expected 
to be an analytic function of T at any value of the tem- 
perature, and it can vanish only at T = 0. There is 
no Kauzmann temperature here or in any similar generic 
model. 



C. The notion of Active temperature 

Notwithstanding all of the above, the system under 
study can slow down so much that upon reducing the 
temperature one has to wait for a very long time before 
equilibrium is reached. When the relaxation times are 
already very large, say at an initial temperature Tj, any 
rapid decrease in the temperature of the heat bath to a 
final temperature Tf may result in a very lethargic re- 
sponse of the system, which may keep the concentrations 
of various quasi-specics at values which are consistent 
with T; rather than Tf . It is then perfectly legitimate to 
introduce the notion of a Active temperature Tg c ~ T, 
as long as one is satisfied with short observation times. 
For longer and longer times the system will exhibit the 
process of aging, and in particular Tfi c will converge to 
Tf , reaching there with certainty if given enough time. 

IV. SLOWING DOWN AND ENTROPIC 
SQUEEZE 

The aim of this section is to explain the most impor- 
tant aspect of glass formation, i.e., the extreme slowing 
down in relaxation to equilibrium when the temperature 
is lowered. The riddle is as follows: the natural time 
scale is determined by the molecular jitter due to ther- 
mal motion. This time is typically of the order or 10 _12 s 
at room temperature. Glassy dynamics exhibits relax- 
ation time T e of the order of seconds, or hours, sometimes 
years. How is it that such a huge gap in time scale is ob- 
tained without geometrical obstruction (as is the main 
theme of this paper)? We will explain that the issue is 
cntropic squeeze, or the failure of entropy or 'the num- 
ber of available paths' to overcome the necessary energy 
climb required for relaxation. 

In 0, [32j it was argued that the relaxation time can 
be predicted if one knew the typical scale £ that sepa- 
rates 'liquid-like' quasi-species. In other words, having 
the concentration Ci of large particles in pentagons and 
small ones in heptagons one introduces a typical scale £ 
by 

since the system has two dimensions. To connect be- 
tween the relaxation time r e and the length scale £ it 
was asserted that for the viscous fluid there exists a free 




FIG. 9: (Color online). An example of a phase obtained by 
simulation at low temperatures when hexagons appear in the 
phase of Fig. [8] Note the tendency of hexagons to clump 
together, a tendency that is ignored in our approximate esti- 
mate of the configurational entropy. 

energy of activation AG*(T) associated with the relax- 
ation event, 

T c =T exp(AG*(T)/T) , (11) 

where To is a microscopic time scale of the order of a 
single particle vibration time. The free energy of activa- 
tion is estimated as the number of Voronoi cells N* (T) 
involved in the relaxation event, times the (temperature 
independent) chemical potential per cell A/x, AG*(T) w 
N*(T)A/i. The number N* depends on whether the re- 
laxation event is a 1-dimensional [3^ | or 2-dimensional 
event 0. In the first case N*(T) w £(T)_/vT2 while in 
the second N*(T) « 7r£ 2 (T)/4ft, where Q is the mean 
area of a Voronoi cell. We end up with the predictions 

r e = t exp(^(T)A/i/V^T) ,1-d event , (12) 
T e = t cxp« 2 (T)A^/4fjT) ,2-d event . (13) 

These predictions were shown to fit the simulation data 
very well [H, |32| . Here we want to explain the fundamen- 
tal origin of these formulae. 

The glass transition and the associated slowing down 
take place in the range of temperatures around T± where 
the liquid-like quasi-species deplete rapidly. It is advanta- 
geous theoretically to focus on the range of temperatures 
around T2 where the hexagons deplete quickly, cf. Fig. [91 
since then we have a smaller number of quasi-species to 
take into account, while the fundamental phenomenon of 
cntropic squeeze is not very different. So think about a 
situation when the majority of quasi-species are pairs of 
glass-like pentagons and heptagons, and set up the en- 
ergy scale such that these pairs have energy zero. Next 
consider a temperature T where the equilibrium concen- 
tration c of hexagons is small, c< 1. Set up the energy 
units such that each pair of hexagons (one with small par- 
ticle and one with large particle) has an energy A = 1. 
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Accordingly the energy of this configuration is 

E = cN . (14) 

One the other hand, the entropy of this configuration can 
be estimated as the logarithm of the number of ways the 
hexagons can be distributed, which is 

S ps -NcAtlc = Ncbx(l/c) . (15) 

As a result we can compute 

— — = — = — lnc— 1m— Inc. (16) 

dE N dc v ' 

Remembering the thermodynamic identity dE/dS = T 
we then conclude that the concentration of hexagons sat- 
isfies the relation 



c = exp(-l/T) . 



(17) 



Accordingly we conclude that in two-dimensions the av- 
erage distance £ between hexagons, which is £ ps l/y/c, 
satisfies 



i « exp(l/2T) 



(18) 



Whenever a pair of hexagons is created, with high 
probability, this will be undone sometime in the future 
(by the inverse operation). Once they are separated in 
space by the distance £, we need a number of flips which 
is at least of the order of £, but maybe many more, in or- 
der to annihilate a pair. This is the fundamental process 
of relaxation at temperature T which we now proceed to 
estimate. In other words, we estimate how many flips 
are typically needed in order to get rid of one pair of 
hexagons, and we will show that the answer is super- 
Arrhcnius. 

We first note that because the concentration c is small, 
any random flip will create a pair of hexagons and in- 
crease the energy by unity. Only flips that eliminate a 
pair of hexagons (to create a heptagon and pentagon) re- 
duce the energy by unity. We denote the probability of 
such a rare event by 7. It will be crucial that 7 depends 
on the temperature as we will explain below. Out of all 
the other flips, assume that a fraction a of flips does not 
change the energy (it has to involve a hexagon and one of 
a pentagon-heptagon pair). What is then the best way to 
move one hexagon a distance £ until it can annihilate its 
counterpart? Assume that the path takes m flips where 
the energy increases by one unit (a typical flip), is energy 
neutral over k flips, and then goes down in energy in m 
steps, with the constraint 



2m + k > £ . 

Of course, these events can take place in any order. 
Clearly, we have a competition between the number of 
ways to arrange such a path and the energy barrier that 
needs to be surmounted, i.e., we need to sum up over all 
m and k the expression 



P(m, k) 



(2m)! 
m! m! fc! 



Note that all the factors are smaller than 1 and therefore 
the sum is well approximated by the largest term, which 
occurs for m = £/2 and k = 0. We thus get 



P= P(m,k) ~ [(1 - a -j)je 

2m+fc>£ 



-l/Tif/2 



(21) 



Note now that a is just some number smaller than 1, 
but more importantly, 7, which is the probability to find 
an energy-lowering move is proportional to the density 
of defects. Indeed, if any flip is done at a link in the 
triangulation in which neither of its ends is a defect the 
energy will go up. (In fact, this reasoning also shows that 
a goes to zero when the temperature goes to 0.) 

Now, the density of defects is c = exp(— 1/T) = 7 so 
that we get finally, 



P~ (e~ 



l/T e -l/T)exp(l/2T)/2 



The time we need to wait to see the event is proportional 
to the inverse of this probability, or 

r e ex exp (£/T) (22) 

which is precisely Eq. (|12|) . This is the one appropri- 
ate for 1-dimcnsional relaxation events as assumed here 
and the generalization to other dimensions is obtained by 
modifying the relation between the density and £. The 
non-Arrhenius nature of the relaxation time is due to the 
strong temperature dependence of £(T) (cf. [H]), which 
in turn is due to the fast disappearance of a class of quasi- 
species. This reduction of the number of quasi-species is 
responsible for the entropy squeeze. 



(1 - a - 7 ) m a fe 7 m e 



V. SUMMARY AND CONCLUSIONS 

Glass forming systems with hard cores can get jammed 
because of geometric constraints. In this paper we ar- 
gued, on the basis of a generic example, that when the 
potentials are soft, the spectacular slowing down associ- 
ated with the glass transition is in a sense more interest- 
ing, since it does not occur due to geometric jamming. 
Such system never lose their ergodicity, and any config- 
uration can be reached from any other in a polynomial 
(in N) number of steps, even though the number of con- 
figurations is exponential in the number of particles. We 
demonstrated explicitly that the configurational entropy 
in such systems if finite at any temperature, and thus nei- 
ther the Kauzmann temperature nor the Vogel-Fulchcr 
(19) formula can be taken seriously. Both are the result of 
an extrapolation which is not fundamental. Finally we 
addressed the question of what is the reason for slowing 
down, and what is the mechanism of its super-Arrhcnius 
temperature dependence. We showed by an explicit cal- 
culation for a generic relaxation step that near the glass 
transition when the concentration of some quasi-species 
becomes very small, the entropic squeeze results in the 

(20) 

inability of the entropic counting of paths to balance the 
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energetic barriers, leading therefore to relaxation times 
that depend on the temperature much faster than ex- 
pected from the Arrhcnius form. 

In summary, we focused on the topological properties 
of generic glass forming systems, to clarify some funda- 
mental issues which are not always clear in the litera- 
ture. Needless to say, much of the interest in glass form- 
ing systems, including their mechanical properties, calls 
for understanding further issues, including metric issues 
that are outside the scope of this paper. For some recent 



thoughts on these subjects we refer the reader to [l9l.l3ll]. 
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